Simulations of grafted polymers in a good solvent 
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We present improved simulations of three-dimensional self avoiding walks with one end attached 
to an impenetrable surface on the simple cubic lattice. This surface can either be a-thermal, having 
thus only an entropic effect, or attractive. In the latter case we concentrate on the adsorption 
transition, We find clear evidence for the cross-over exponent to be smaller than 1/2, in contrast 
to all previous simulations but in agreement with a re-summed field theoretic e-expansion. Since 
we use the pruned-enriched Rosenbluth method (PERM) which allows very precise estimates of the 
partition sum itself, we also obtain improved estimates for all entropic critical exponents. 



I. INTRODUCTION 

Single polymers grafted to a plane impenetrable wall 
have been studied intensely for a long time. If the surface 
is energetically neutral, the main effect is a change in the 
critical exponent governing the scaling of the partition 
sum. For an attractive surface the situation is more in- 
teresting, as there is a second order adsorption transition 
at a finite strength of the attraction 

Of particular interest is the cross-over exponent at this 
transition. Early simulations 1, 2] had given <fr > 0.5, but 
it was claimed in Q that this was due to finite size correc- 
tions to scaling, and that <fi is actually very close to 1/2 
(the best estimated being <j> = 0.496 ± 0.004). Basically 
the same conclusion was reached recently in |4| , where an 
even smaller error bar was obtained, <fi = 0.5005 ±0.0036, 
suggesting that 0=1/2 exactly. Since 0=1/2 also 
in d = 2 and d > 4, this would mean that <f> is supe- 
runiversal, as it is for branched polymers [5J, |6J. But 
a completely different picture was drawn in another re- 
cent paper pj, where it was claimed that <f> — 0.59. To 
add to the confusion, we should cite field theoretic re- 
sults. The e-expansion with e = 4 — d predicts @, 
(f> = l/2-e/16+[167r 2 -39]e 2 /512 + .... This would give 
(f) = 0.68, if higher order terms were simply omitted, but 
Pade-Borel summation gives 6 = 0.483 101. On the other 
hand, fixed dimension (massive theory) renormalization 
group calculations [Hj,[ll( give 4> — 0.52. All these calcu- 
lations have errors which are hard to pin down, and the 
authors of [ljj preferred the value 0.52 over 0.483. Any- 
how, these renormalization group calculations strongly 
suggest that <f> is not superuniversal. 

It is the purpose of the present work to clarify the sit- 
uation by means of much more precise simulations. We 
will find that indeed (j> is definitely smaller than 1/2 (as 
predicted by the resummed e-expansion, but not by the 
supposedly more reliable massive field theory approach) . 
In addition, wc shall also provide precise estimates for 
the location of the adsorption transition and for the en- 
tropic exponents. The latter will be done both for ther- 
mal surfaces (at the adsorption transition point) and for 
a-thermal surfaces. 

We model the polymer by a self avoiding walk (SAW) 
of N — 1 steps on a simple cubic lattice with restriction 



z > 0. There is an energy — e for each monomer (site) 
at z = 0, the first monomer is located at x = 0. The 
temperature is taken to be T = 1 , so that the Boltzmann 
factor for each contact with the surface is q = e e . The ad- 
sorption transition is at q = q c > 1. For the simulations 
we use the pruned-enriched Rosenbluth method (PERM) 
This is a recursively (depth-first) implemented 
sequential sampling algorithm with re-sampling [l4| . It 
is similar to the algorithm used in ||, but it is faster 
and much easier to use. Its main advantage over con- 
ventional Markov chain Monte Carlo (MC) methods is 
that it gives very precise estimates of the partition sum, 
without any need for thermodynamic integration or the 
like. To minimize statistical errors and speed up the al- 



ieea up tj 

gorithm, we use Markovian anticipation [Til lli| . We 
simulated w 6.6 x 10 8 walks with N = 8000 for q = 1 
and « 9.5 x 10 8 walks at q « q c . In both cases, about 1.6 
% of these walks were strictly independent. Altogether 
this needed w 5000 hours CPU time on fast (3GHz) PCs. 
During the runs with q > 1, results at slightly different 
values of q were obtained by re-weighing on the fly, so 
that one run made with q w q c gave results at three close 
values of q. The critical point q c was then found by in- 
terpolation. 

The partition sum is written as Z\(N,q) — 
CNmQ m , where CV m is the number of configurations 
with to contacts with the wall, Zi k = for k = 1, ... to, 
and the subscript '1' indicates that one end is grafted. 
For q < q c it scales as 



%i q) ~ n N 1 



(1) 



with [i and 71 independent of q, but with a g-dependent 
prefactor. 

Near the adsorption transition, Z\{N,q) should scale 
as 1 1 



Z,(N, q) ~ ^N^-^Kq - q c )N% 



(2) 



where ty(z) is analytic for finite z and lim z ^_ 00 ^f(z) is 
finite and positive. Notice that the dominant exponen- 
tial growth of Z\ (JV, q) with N is the same as for ordi- 
nary SAWs, as long as q < q c . Taking the derivative of 
In Z\ (N, q) with respect to q and setting q = q c there- 
after, we obtain for the average energy exactly at the 
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critical point 



E N (q c ) = (em) ~ JV* 



(3) 



Taking two derivatives we would obtain a scaling ansatz 
for the specific heat which is often used to estimate <f> 
and other aspects of the critical behavior. We will not 
use it in the present paper, since Eqs.© and (0) give 
much more precise results, as already found in [Tj- Ljl - We 
refer to |(J (which deals with the analogous problem for 
branched polymers) for a detailed discussion. Notice that 
Eq.© cannot be used with Markov chain MC methods, 
since the latter do not give simple and precise estimates 
of Zx(N,q) itself. 



II. RESULTS: A-THERMAL WALLS 

From previous simulations 0, Il7| we know that /i w 
4.68404, but the present simulations have higher statis- 
tics, therefore we should first estimate the critical expo- 
nent 71 in such a way that neither fj, nor the unknown 
prefactor in Eq.QJ affect the value. For this we form the 
triple ratios 



.... 41nZi(JV)-31nZi(JV/3)-lnZi(3^) 
7i,e ff (iV) = 1 + — 

(4) 

(with Zi(N) = Zi(N,l)), which should tend to 71 as 
N — > 00. The leading corrections to Eq.QJ (and thus also 
to Eq.Q) should scale as 1/A A with the same exponent 
A sa 0.5 holding also for SAWs in absence of a surface 0, 
ITsl IT^ | . This would suggest that we should get a straight 
line when plotting 71 (N) against 1/y/N. Unfortunately 
this is not true, due to the presence of very large analytic 
corrections cx 1/N, 1/N 2 , .... As seen from Fig. 1, they 
shift the effective correction to scaling exponent to A D g w 
0.7, and the extrapolation to N — > 00 gives 



71 = 0.6786 ±0.0012 



(5) 
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FIG. 1: Effective exponents defined by Eq.Q, plotted against 
1/iV ' 7 . 
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FIG. 2: The combination lnZi(A) - iVln/i + (1 - 71) ln/V, 
with 71 determined from Fig. 1 and for three candidate values 
of (i, plotted against 1/N ' 7 . Error bars are plotted only for 
the central curve, but errors are the same (within the precision 
of drawing) for all three curves. 



This is in good agreement with the best previous MC 
estimate 0.679 ± 0.002 and with the field theoretic 
result 0.680 pjj. 

Using this value of 71 , we next show in Fig. 2 the quan- 
tity lnZi(iV) -aN+ (1 - 71) In N for different values of 
a, plotted again versus l/iV A . Again these data should 
extrapolate to 1/N — > with a finite slope, if a = ln/i. 
From Fig. 2 we see that 

l/fi = 0.21349098±0.00000005 (ji = 4.6840386±0.0000011) 

(6) 

where the error includes the uncertainty in the estimate 
of 71. This is the most precise value of the critical fu- 
gacity of SAWs on the simple cubic lattice published up 
to now. It is about 80 times more precise than the best 
estimate from exact enumerations |2l). 

In addition to the partition sum for SAWs grafted at 
one end to the surface, we can also study the partition 



sum Z\\{N) where both ends are constrained to have 
z — 0. In analogy to Eq.JTJ, it should scale as 



Z n (N,q= 1) - ^jV 711 ^ 1 



(7) 



The exponent 711 should be related to previous expo- 
nents by the Barber relation [2^ 



7 - 271 + 711 + 







(8) 



where 7 and v arc the entropy and Flory exponents (i? ~ 
N v ) for ungrafted SAWs in the bulk. Values of Z n (N) 
are obtained simply by summing over those walks for 
which the A-th monomer has z = 0. Using the above 
estimate of fi and different candidate values of 711, we 
plot in Fig. 3 the analogous quantity to that shown in 
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FIG. 3: Analogous to Fig. but for SAWs grafted at both 
ends to the surface, and with the candidate values of fi re- 
placed by candidate values of 711. In order to reduce statisti- 
cal errors, the data are binned with bin width AN/N ~ 0.01. 
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FIG. 4: (x 2 N + y 2 N )/N 2 " (upper curve) and {z 2 N )/N 2u (lower 
curve) against 1/yN . Error bars are smaller than the thick- 
ness of the lines. 



Fig. 2 for singly grafted chains. We find 
711 = -0.390 ±0.002 , 



(9) 



to be compared to the previous MC estimate —0.383 ± 
0.005 3] and to the field theoretic prediction -0.388 [l0| . 
Using the most precise previous estimates of 7 and v from 
H3 (7 = 1.1573 ±0.0002 and v = 0.58765 ± 0.00020), we 
see that the Barber relation is indeed perfectly satisfied, 



7 - 271 + 711 + 



-0.0023 ±0.0031. 



(10) 



Root mean square end-to-end distances should scale 
with the same Flory exponent as in the bulk, but with 
different prefactors and with different corrections to scal- 
ing. We measured both the components parallel and per- 
pendicular to the wall for singly grafted SAWs. This time 
the corrections to scaling were ~ 1/VlV as expected, i.e. 
there are much smaller analytic corrections. Results are 
shown in Fig. 4, where we divided averaged square dis- 
tances by N 2u and plotted them against l/y^N. The 
ratio (z 2 ^) / (x 2 N + y 2 N ) increases with N as found also 
in [23[, but it converges for TV — > 00 to a finite value, 
0.938 ± 0.002. In contrast to |2^| we see no indication 
that either (z%) or (x 2 N + yjy) scales with an exponent 
different from the bulk Flory exponent. 



III. RESULTS: ATTRACTIVE WALLS 

According to Eq. J2J , the exponential growth of the par- 
tition sum with N is the same at q c as at q = 1, i.e. we 
do not need the triple ratio Eq.(0J to estimate 7* (.AT). 
Instead we can use 



7? off (AU) = l + 



\n[Z 1 (2N,q)/Z 1 (N/2,q)/ f i 3N / 2 } 
In 4 



(11) 
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FIG. 5: Effective exponents 7 l off (A r , g), plotted versus 1/yN, 
for several values of q. 



where the ratio between the two partition sums elimi- 
nates the unknown prefactor. The critical point is char- 
acterized by the fact that eS (N, q) diverges for q > q c 
when N — > 00 , converges slowly to 71 for q < q c , and con- 
verges to a constant value larger than 71 exactly at q c . 
Values of 7* off (iV, q) obtained by means of EcL lfTO are 
shown in Fig. 5. As in the next figures to follow, we plot- 
ted it against 1/yN since there was no different unique 
value of A which fitted all observables, and A = 1/2 was 
not worse overall than other values. 

From Fig. 5 we see that q c « 1.3307 and 75 s ~ 1.23, in 
good agreement with the estimates of 3J. More precise 
values will result by combining also the information from 
other observables. 

The next quantity we looked at is the average en- 
ergy, or rather the average number of sites in contact 
with the wall. Plotting simply (m)/N^ would suggest 
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FIG. 6: Average number of contacts with the wall, divided by 
■\/N. This would suggest <f> ~ 1/2, if we would assume that 
corrections to scaling are negligible for N > 4000. 
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FIG. 8: Difference between the effective entropic exponent for 
polymers grafted at one end, and the exponent for polymers 
grafted at both ends, again plotted versus 1/yN. 
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FIG. 7: Effective cross-over exponent as defined in the text, 
plotted versus 1/viV. 
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FIG. 9: Ratio (z 2 )/(x 2 +y 2 } between mean square end-to-end 
distances perpendicular and parallel to the wall. 



4> = 1/2 and q c = 1.331, see Fig. 6. But this does 
not take into account the fact that corrections to scal- 
ing are very large and should still be important even for 
the largest N. Indeed, when defining an effective ex- 
ponent by (t) cS (N,q) = (ln4) _1 \n[E 2N (q) /E N/2 (q)] and 
plotting it against l/\f~N (Fig. 7), we see that it extrap- 
olates clearly to a value < 1/2, and that q c is closer to 
the value 1.3307 found from Fig. 5. 

Instead of the effective exponent 7^ cff (iV, q) defined in 
complete analogy to Eq. lfTT|) . we show in Fig. 8 the dif- 
ference cS (N, q) — tJ x eS (N, q). This is obtained from 
the ratio between the partition sums for singly and dou- 
bly grafted polymers, and is independent of the precise 
value of /z. Finally, we show in Fig. 9 the ratio between 
the mean square end-to-end distances perpendicular and 
parallel to the wall. 

Comparing Figs. 5 to 9, we see that all of them (except 



of course Fig. 6) give the same value for q c . Thus we 
obtain our precise estimate 

q c = 1.33065 ±0.00010. (12) 

This is in agreement with the value 1.3310 ± 0.0003 of 
@, and with the recent estimate 1.334 ± 0.026 of 0. For 
the anisotropy of infinitely long chains at the adsorption 
transition this gives (z 2 )/{x 2 + y 2 ) = 0.3845 ± 0.0010. 
The best estimates of the critical exponents are then 

4> = 0.484 ±0.002, 75 s = 1.226 ± 0.002, (13) 

and 

7l s - 7 f 1 = 0.519 ± 0.003. (14) 

These values are in less agreement with previous esti- 
mates. First of all, the entropic exponents are slightly 
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smaller (by about 2 standard deviations) than the values 
given in [3j, although these were already lower than all 
previous MC estimates. While 7* agrees roughly with 
the field theoretic prediction 1.207 of |10|. 7^ is quite a 
bit larger (0.707 as opposed to 0.666). Secondly, and 
more importantly, the cross-over exponent is now clearly 
less than 1/2, by some 8 standard deviations. The clos- 
est previous MC estimate was 0.496 ± 0.004 |3j, which is 
off by three sigma. The estimate 0.5005 ± 0.0036 of 4) 
seems clearly excluded by our data. It is based on very 
high statistics of very short (N < 200) chains, which 
suggests that the corrections to scaling were not taken 
fully into account in Q. As mentioned above, (j) < 1/2 
is predicted by the first order term of the e-expansion, 
but not when terms up to 0(e 2 ) are included and not 
by fixed dimension renormalization group methods. But 
the resummation of the e-expansion done in [Toj gave 
<j> = 0.483, in surprisingly good (and presumably fortu- 
itous) agreement with our result. At least, our estimate is 
comfortably larger than the first order epsilon expansion 
result, 0> l/2-e/16 = 0.4375. 



IV. CONCLUSION 



In this paper we have presented Monte Carlo simula- 
tions of single grafted 3-d polymers in a good solvent, 
both for attractive and for a-thermal walls. The sample 
seems to be the biggest studied so far, both as concerns 
the number of chains simulated and their lengths. By 
using PERM which gives precise estimates of the parti- 
tion sum, we could use the partition sum itself (instead of 
the specific heat) to locate the critical adsorption point 
and the critical exponents. Our estimates continue the 
decrease with increased statistics observed already in Q • 
Our most interesting result is that the cross-over expo- 
nent is clearly less than 1/2, in contrast to all previous 
simulations and to the best estimates from field theory 
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prior to publication. 
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